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It was shown that in the canonical ensemble the simple exactly soluble statistical model of nuclei 
decay into nucleons, which is a limiting case of the statistical multifragmentation model, predicts 
the nuclear first order phase transition associated with the Helmholtz free energy different from the 
first order phase transition of the liquid-gas type associated with the Gibbs free energy. The main 
thermodynamic properties of this phase transition were explored on the basis of the method of the 
thermodynamic potential and its first and second derivatives in the thermodynamic limit. It was 
established that the thermodynamic potential F is a piecewise smooth function and its first order 
partial derivatives with respect to variables of state are piecewise continuous functions. At the points 
of phase transition, the energy in the caloric curve is discontinuous at the constant temperature and 
fixed values of the specific volume, while the pressure and the chemical potential in the equations 
of state are discontinuous at the constant specific volume and fixed values of the temperature. 
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I. INTRODUCTION 



The heavy-ion collision experiments at RHIC, the LHC and the future FAIR and NICA projects are designed 
to establish the QCD phase diagram. The first order QCD phase transition, which separates the hadronic 
^ 1 phase from the deconfined quark-gluon phase below the critical endpoint, is thought to be represented by a 
coexistence line on the phase diagram in the plane of temperature T and baryon chemical potential fj, 
Such a first order phase transition, which is represented by a line in the T — fi plane, is a liquid-gas type one 
associated with the Gibbs free energy, which is a thermodynamic potential of the isobaric ensemble [8|. It 
■ is accepted that the nuclear multifragmentation, i.e. the disintegration of nuclei into fragments observed in 
(*C) [ intermediate-energy nuclear reactions, is also the first order phase transition of the nuclear liquid-gas type (see, 
lO ■ e.g., ref. Q, and the references cited therein). However, on the theoretical side there is an ambiguity about 
C ' ' the existence of discontinuity of energy (plateau) in the nuclear multifragmentation caloric curve in different 
statistical ensembles. One class of the statistical multifragmentation models (SMM) predicts a plateau in the 
\Q ■ caloric curves in the canonical and microcanonical ensembles [Tol - [T3j . but another class of SMM predicts a 
plateau in the caloric curves in the isobaric ensemble [Til [l5j , which is a typical signature of the nuclear liquid- 
^1 ■ gas phase transition. Another example of the first order phase transition of the liquid-gas type which usually 
' m 1 1 happens in nuclear matter is the phase transition taking place in the interacting hadronic nuclear matter at low 
> ■ temperatures fj. fl6l [l7j . 

The first order phase transition of the liquid-gas type associated with the Gibbs free energy is defined by the 
cusp of the Gibbs free energy per particle (the chemical potential) and the corresponding jump discontinuities 
of its first derivatives, i.e. the specific volume and the entropy per particle, at constant values of pressure and 
temperature. Namely, the Gibbs free energy per particle is a piecewise smooth function and its first order partial 
derivatives with respect to variables of state of the isobaric ensemble (T,p), i.e., the entropy per particle and 
the specific volume are the piecewise continuous functions. For the liquid-gas phase transition the energy in the 
caloric curve is discontinuous in the isobaric and the grand canonical ensembles at fixed values of the pressure 
and the chemical potential, respectively, and it is continuous in the canonical and microcanonical ensembles at 
fixed values of the specific volume. However, the specific volume in the isotherms is discontinuous in the isobaric 
and the canonical ensembles at fixed values of the temperature. The phase diagrams for the liquid-gas phase 
transition are represented by the coexistence lines in the planes T — fj, and T — p and the coexistence areas in 
the plane of temperature and specific volume (density of particles) and in the plane of temperature and entropy 
per particle. 

In the Ehrenfest classification scheme the first order phase transition is defined as a phase transition which 
can be associated only with the Gibbs free energy of the isobaric ensemble [H, [l9| . Such a phase transition is a 
liquid-gas type one. In refs. [I(|[2lj|, the Ehrenfest definition of the first order phase transition was generalized to 
any thermodynamic potential. In this case, it is the one associated with a finite discontinuity in one or more first 
derivatives of the appropriate thermodynamic potential with respect to its variables of state. Some models which 
predict the first order phase transition in the canonical and microcanonical ensembles, the thermodynamical 
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potentials of which arc different from the Gibbs free energy, can be seen in refs. [22h26||. 

The present paper proposes a simple statistical model of total nuclei decay into nucleons, which is a limiting 
case of the statistical multifragmcntation model which has the first order phase transition associated with 
the Helmholtz free energy for which the Ehrenfest definition is given in the canonical ensemble. The exact 
analytical results for this model were obtained. They allowed us to find the general thermodynamic properties 
of the first order phase transition associated with the Helmholtz free energy and determine its differences from 
the phase transition of the liquid-gas type. For example, the free energy per nucleon is a continuous function 
at the points of phase transition, but the first order partial derivatives of the thermodynamic potential of the 
canonical ensemble with respect to variables of state, i.e. the entropy per nucleon, the pressure and the chemical 
potential, have jump discontinuities. 

The structure of the paper is as follows. In Section II, we briefly describe basic ingredients of the statistical 
model of total nuclei decay into nucleons. The thermodynamic results for the first order phase transition asso- 
ciated with the Helmholtz free energy and the caloric curve are discussed in Sections III. The main conclusions 
arc summarized in the final section. 



II. THE STATISTICAL MODEL OF THE NUCLEI DECAY INTO NUCLEONS 



A. The general formalism 

The partition function of the statistical multifragmcntation model in the canonical ensemble which describes 
the decay of one nucleus of A nucleons into different nuclear fragments in the volume V at the temperature T 
is actually given by the expression [l3|, [27l - l30j 

( A \ A nk A 

{n k } \k=l I fc=l k ' fe=l 

and 



Uk = 9kVf \2^j 6 ~' (2) 

where Zq = 1, ujk is the partition function of a fragment (nucleus) which has k nucleons, Vf = V — vqA is the 
free volume, Wk is the binding energy of the fragment of k nucleons, is the nucleon mass, v$ = 1/po, po is 
the normal nuclear density and gt is the spin-isospin degeneracy factor. The binding energy of the free nucleons 
Wi = 0. Note that throughout the paper we use the system of natural units, % = c = k B = 1- 

The partition function (fTJ in the limiting case of the SMM, when one nucleus of A nucleons decays into free 
nucleons only, u)% = . . . = ua-i = 0, can be written as 

Let us generalize this limiting case of the SMM to the system of N number of nuclei of A nucleons which 
decay into B free nucleons, B = NA, in a volume V, in contact with a heat reservoir of temperature T. 
Considering the definition ^ and the Maxwell-Boltzmann statistics of nuclei and nucleons, we can write the 
thcrmodynamical potential, the Helmholtz free energy, and the partition function in the canonical ensemble 
(T, V, B) as 



F=-ThxZ B (4) 



and 



Here, the free volume Vf = V — VqB is used. The mean number of the free nucleons (ni) and the mean number 
of the initial free nuclei («a) are given by 

, v B w f / v N 
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The total baryon charge of the system (B) and the mean multiplicities of particles in the system (m) arc 

B = { ni )+A(n A ), (7) 
(m) = {m) + {n A }. (8) 

The first order partial derivatives of the thermodynamic potential F with respect to variables of state of the 
canonical ensemble (T, V, B), i.e. the entropy S, the pressure p and the chemical potential fx can be written as 
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where - 0( z ) is the psi- function. The mean energy of the system is given by 
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The second order partial derivatives of the thermodynamic potential F with respect to variables of state, i.e. 
the heat capacity Cy, the compressibility and the susceptibility Xi can t> e written as 
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The second derivatives allow us to identify the convexity properties of the thermodynamic potential F which 
are related to the thermal, mechanical and particle stability of the system [l9j |. 



B. The thermodynamic limit 



Let us consider the thermodynamic limit: N — > oo, V — > oo and v = V/B = const. Then the free energy per 
nucleon (j4]) reduces to 

-►oo, (16) 
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where w/ = Vf /B = v — vq. The mean number of the free nucleons and the mean number of the free nuclei 
then become 



(m) = b 



i 



l + e 
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where 



a = In — . 



(19) 
(20) 
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The first order partial derivatives of the Hclmholtz potential F with respect to variables of state, i.e. the entropy 
per nucleon ©, the pressure (fTU|) and the chemical potential (fTTj) . in the thermodynamic limit take the form 
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The mean energy per nucleon |T 
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in the thermodynamic limit becomes 
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Finally, the second order partial derivatives of the thermodynamic potential F with respect to variables of state, 
i.e. the heat capacity per nucleon (|13p. the compressibility (|14[) and the susceptibility (|15p. in the thermodynamic 
limit can be written as 
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Note that in this model the thermodynamic limit is taken in accordance with the standard rules of the thermo- 
dynamics and the statistical mechanics (3l| . 



III. THE FIRST ORDER PHASE TRANSITION IN THE CANONICAL ENSEMBLE 



Let us describe the phase transition of the statistical model of the nuclei decay into nucleons in the canonical 
ensemble. This ensemble is characterized by the potential, F(T, V, B) and the physical obscrvablcs correspond 
to its first and second order partial derivatives. The geometry (curvature) of the hypcrsurface F determines the 
physical properties of the system. In the thermodynamic limit at three different values of the variable a, a < 0, 
a = and a > 0, the thermodynamic potential per nucleon (|16[) is a homogeneous function of the first degree 
of the extensive variable of state B: 

F(T,V,B) = Bf(T,v). (28) 

Then in the thermodynamic limit at three different values of the variable a the first and the second order 
partial derivatives of the thermodynamic potential F with respect to variables of state (T, V, B) (fig j) -([2"5 j) can 
be written as 

s = -&(T, v), P = -f' v {T, v), ft = f + pv, (29) 

% = -ft T (T, v), T~~ = C(T, v), — 2 = .C(T, v), (30) 
1 ktv x v 

where f' x {x,y) = (df/dx) y and fx X {x,y) = (d 2 f/dx 2 ) y . These quantities satisfy the differential equation for /, 
the fundamental equation of thermodynamics and the Euler theorem 

df = —sdT — pdv, Tds = de + pdv, Ts = e + pv — [i. (31) 

Moreover, the Legendre transform, / = e — Ts, and the relation for second derivatives, B/x = v/kr, are valid. 
Note that eqs. (|2"8"]) - (f3"Tj) are satisfied separately in the regions a < 0, a = and a > 0. 

The values of a < give the homogeneous phase / of the system which consists of the free nuclei (n\)/B = 
and (n A )/N = 1. The homogeneous phase II of the system is determined by the values of a > in which 
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T(MeV) T(MeV) 
1: (Color online) The thermodynamic potential per nucleon /, the entropy per nucleon s, the pressure p, the 



chemical potential n, the energy per nucleon e, the heat capacity per nucleon c„, the compressibility 1/krv and the 
susceptibility per nucleon B/x as functions of the temperature T at fixed specific volume v for the finite number N of 
nuclei and the thermodynamic limit N = oo. The dotted, dashed and solid curves were calculated by the statistical 
model of nuclei decay into nucleons in the canonical ensemble for the 12 C nuclei with N = 1,3 and N = oo, respectively, 
at the specific volume v = 3vq. The symbols are the phase diagrams in the thermodynamic limit. 



the system is composed by the free nucleons {n\)/B = 1 and (ua)/N = 0. The values of a = in Eq. (l20l) 
determine the line of phase boundary in the T — v space 



v 



Figure Q] presents the behavior of the free energy per nucleon / and its first and second partial derivatives 
with respect to variables of state as functions of the temperature T at fixed specific volume v for the 12 C nuclei 
which decay into nucleons for the finite number N of nuclei and in the thermodynamic limit, N = oo. In the 
calculations we shall assume the typical numerical values of the spin-isospin degeneracy factors adopted in the 
SMM [33], i-e. gi = 4 and qa = 1. Therewith, the binding energy per nucleon for the nucleus of 12 C is taken 
to be Wa/A = 7.680144 McV. The function f(T) is continuous for all T and has a cusp at the point of phase 
transition T = T* where the first order partial derivatives, i.e. the entropy per nucleon s(T), the pressure p(T) 
and the chemical potential //(T), and the energy per nucleon e(T) have jump discontinuities. At this point T* 
the second order partial derivatives, i.e. the heat capacity per nucleon c„(T), the compressibility l/kr(T)v and 
the susceptibility per nucleon B/x{T), have the local extremum (infinite maximum or minimum). Therefore, in 
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conformity with the Ehrenfest definition of the phase transitions the point T = T* is the point of the first order 
phase transition related to the thermodynamic potential F of the canonical ensemble, i.e. the Helmholtz free 
energy. The jump of the entropy per nucleon at the point of phase transition T* is related to the latent heat. 
In this nuclear system two homogeneous phases can be distinguished. In the homogeneous phase (/) at the 
temperature T < T* the system consists of the N free nuclei of 12 C, (tia) = N. However, in the homogeneous 
phase (II) at the temperature T > T* the system is composed of the free nucleons only (m) = B. At T = T* 
the mixed phase is realized. Note that in the canonical ensemble at fixed values of the volume and the number 
of nuclei the pressure of the system is discontinuous at the points of phase transition due to the jump increase 
in the multiplicity of particles which carry the partial pressure ~ T/Vf. The discontinuity of the chemical 
potential is also related to the internal increase in the number of free particles. 




v/v o v/v o 

FIG. 2: (Color online) The thermodynamic potential per nucleon /, the entropy per nucleon s, the pressure p, the 
chemical potential fj,, the energy per nucleon e, the heat capacity per nucleon c v , the compressibility 1/krv and the 
susceptibility per nucleon B/x as functions of the specific volume v at fixed temperature T for both the finite number N 
of nuclei and the thermodynamic limit, N = oo. The dotted, dashed and solid curves were calculated by the statistical 
model of nuclei decay into nucleons in the canonical ensemble for the 12 C nuclei with N = 1, 3 and N = oo, respectively, 
at the temperature T = 18.4 MeV. The symbols are the phase diagrams in the thermodynamic limit. 

Figure [2] presents the behavior of the free energy per nucleon / and its first and second partial derivatives 
with respect to variables of state as functions of the specific volume v at fixed temperature T for the 12 C nuclei 
which decay into nucleons for the finite number N of nuclei and in the thermodynamic limit, N = oo. In the 
thermodynamic limit the function f{v) is a piecewise smooth function on a closed interval / around the point of 
phase transition v — v* , i.e., it is a one- valued continuous broken-line function with the point of discontinuity of 
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FIG. 3: (Color online) The Helmholtz free energy per nucleon f(T, v) as a function of the temperature T and the specific 
volume v for the statistical model of nuclei decay into nucleons in the canonical ensemble in the thermodynamic limit 
for the 12 C nuclei. The solid lines are the line of phase transition on the surface f(T, v) and its projection onto the plane 
T — v, i.e. the phase diagram. 

the first derivative at v = v* . The graph has a cusp and there is no derivative. The first order partial derivatives 
of the thermodynamic potential per nucleon with respect to variables of state, the entropy per nucleon s(v), 
the pressure p(v) and the chemical potential p(v) are single-valued piecewise continuous functions for all v. At 
the point v — v* the first order partial derivatives, i.e. the entropy per nucleon s(v), the pressure p(v) and the 
chemical potential p(v), and the energy per nucleon e(v) have jump discontinuities. However, the second order 
partial derivatives of the thermodynamic potential per nucleon with respect to variables of state, i.e. the heat 
capacity per nucleon c v (v), the compressibility l/kT{v)v and the susceptibility per nucleon B/x(v), have the 
local extremum (infinite maximum or minimum) . Therefore, the point of discontinuity of the first derivatives at 
v = v* is the point of the first order phase transition related to the thermodynamic potential F of the canonical 
ensemble, i.e. the Helmholtz free energy. The jump of the entropy per nucleon at the point of phase transition 
v* is related to the latent heat. In the nuclear system described by the statistical model of the nuclei decay into 
nucleons in the canonical ensemble two homogeneous phases are realized. In the homogeneous phase (/) at the 
specific volume v < v* the system consists of the N free nuclei, (ua) = N. However, in the homogeneous phase 
(II) at the specific volume v > v* the system is composed of the free nucleons, (ni) — B. At v — v* the mixed 
phase is realized. 

Figure [3] shows a picture of the surface of the thermodynamic potential per nucleon / as a function of the 
variables of state (T,v) near the line of phase transition which is determined by eq. ([52]) . This line is a line of 
fracture of the surface / as in all its points the function f(T, v) is continuous but undifferentiated. It defines 
the first order phase transition associated with the Helmholtz free energy, and its projection onto the T — v 
plane gives the corresponding phase diagram. 

More concretely, the phase diagrams for the first order phase transition of the statistical model of the nuclei 
decay into nucleons in the canonical ensemble (T, v) in the thermodynamic limit are depicted in Fig. |4j The 
phase diagram T—v for the first order phase transition in the canonical ensemble is represented by the continuous 
coexistence curve which begins at the critical point (T c = oo, v c = vq) and finishes at the point (T = 0, v = oo). 
The coexistence curve T — p crosses the p axes at the point (T = 0, p = 0) and finishes at the critical point 
(T c = oc, p c — po). Along these coexistence curves, the functions s(T, v), p(T, v) and p(T, v), are discontinuous. 
Therefore, in the phase diagrams T — s, T — p, T — p,, v — s, v — p and v — p, we have the coexistence regions 
instead of the coexistence lines because the variables s, p and p are undefined at the points of phase transition 
(T*,v*). We have two homogeneous phases in the system: the phase (I) of the free nuclei and the phase (//) of 
the free nucleons. The mixed phase is defined by the coexistence lines on the phase diagrams T — v, T — p and 
the coexistence areas on the phase diagrams T — s, T — p and T — p. In contrast to this phase transition, the 
first order phase transition of the liquid-gas type associated with the Gibbs free energy G is described by the 
coexistence lines on the phase diagrams T — p and T — p and the coexistence areas on the phase diagrams T — v, 
T — p Q. From the phase diagrams shown in Fig. [3l the Gibbs phase rule for the first order phase transition 
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FIG. 4: (Color online) The phase diagrams T — v, T — p (top panels), T — a, T — p, T — p (mean panels) and v — s, 
v —p, v — p (bottom panels) for the first order phase transition of the statistical model of the nuclei decay into nucleons 
for the s 2 C nuclei in the canonical ensemble. Roman numerals denote the homogeneous phases (7) and (77), the shaded 
areas correspond to the mixed phase and the lines are the coexistence curves. 

associated with the Helmholtz free energy F of the canonical ensemble can be defined by 

Tj = T II} vj = v n , (33) 

where the roman numerals (7) and (77) denote the homogeneous phases in equilibrium. It is characterized by 
the equality of temperature and specific volume across phase boundaries. In contrast to this, the Gibbs phase 
rule for the first order phase transition of the liquid-gas type associated with the Gibbs free energy is represented 
by the equality of temperature, pressure and chemical potential across phase boundaries [33j . 

Figure [5] presents the caloric curve, the temperature T as a function of the multiplicity, the equation of state 
or the isotherms, i.e. the dependence p — p at fixed temperature T, and the dependence p — p in the canonical 
ensemble for the first order phase transition of the statistical model of the nuclei decay into nucleons. The 
excitation energy per nucleon is defined by the equation, e* — e(T) — e(0), where e(0) = —Wa/A is the energy 
per nucleon in the ground state. At the point of phase transition at temperature T = T* the excitation energy 
per nucleon e* in the caloric curve has a jump discontinuity in the canonical ensemble at fixed specific volume 
v or baryon density p. The curve e*(T) is a piecewise continuous function. The multiplicity of particles in 
Fig. [SJb) has a jump discontinuity at the point of phase transition. In the phase (I) it is equal to (m) = N, but 
in the phase (II) the multiplicity is (m) = B. The pressure p and the chemical potential p in the isotherms and 
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FIG. 5: (Color online) (a) The temperature T as a function of the excitation energy per nucleon e* (the caloric curve) 
and (b) the temperature T as a function of the multiplicity (m) for the statistical model of the nuclei decay into nucleons 
for the 6 2 C nuclei in the canonical ensemble for V = 3 and in the thermodynamic limit at the specific volume v = 2.5wo- 
(c) The isotherms and (d) the dependence of the chemical potential /j, on the baryon density p in the canonical ensemble 
for TV = 3 and in the thermodynamic limit at the temperature T — 21.4 MeV. 



the equation p — p, respectively, are discontinuous in the canonical ensemble at the points of phase transition. 
See Fig. |5jc) and (d). Note that in the case of the first order phase transition of the liquid-gas type in the 
canonical ensemble, the excitation energy per nucleon e* in the caloric curve is continuous at fixed values of the 
specific volume and the specific volume v in the isotherms is discontinuous at constant values of the pressure 
p and the chemical potential p [8[. In the canonical ensemble the liquid-gas phase transition is defined by the 
plateaus in the isotherms and the equation p — p at fixed temperature T. 

Summarizing, we have found that the phase transition of the statistical model of the nuclei decay into nucleons 
in the canonical ensemble is the first order phase transition defined by the Helmholtz free energy, instead of 
the Gibbs free energy which defines the first order phase transition of the liquid-gas type. The first order 
phase transition associated with the Helmholtz free energy is defined by the piecewise smooth function of the 
Helmholtz potential per nucleon / and the piecewise continuous functions of the first order partial derivatives 
of F with respect to variables of state, the entropy per nucleon s, the pressure p and the chemical potential p. 
At the points of phase transition given by the variables of state (T, v) the potential / is a continuous function 
which has a cusp, however, the first derivatives s, p and p (order parameters) have jump discontinuities and 
the second order partial derivatives c v , 1/hrv and B/x have an infinite maximum or minimum. The excitation 
energy per nucleon e* in the caloric curve is discontinuous in the canonical ensemble. 



IV. CONCLUSIONS 



To conclude, in this paper the statistical model of the nuclei decay into nucleons in the canonical ensemble 
as the limiting case of the statistical multifragmentation model was introduced. The exact analytical results 
for both the finite system of N nuclei and the infinite system in the thermodynamic limit were obtained. It 
was revealed that the statistical model of the nuclei decay into nucleons in the canonical ensemble in the 
thermodynamic limit has the first order phase transition associated with the Helmholtz free energy instead of 
the Gibbs free energy. This phase transition was investigated on the basis of the method of the thermodynamic 
potentials and their first and second derivatives. The main thermodynamic properties of this first order phase 
transition in the canonical ensemble were found. These properties are as follows. The Helmholtz free energy 
per nucleon is the piecewise smooth function and its first order partial derivatives with respect to variables 
of state, i.e., the entropy per nucleon, the pressure and the chemical potential, are the piecewise continuous 
functions. At the points of phase transition, the free energy per nucleon is a continuous function which has 
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a cusp both as a function of v at fixed T and as a function of T at fixed v; however, the first order partial 
derivatives of the free energy with respect to variables of state, i.e. the entropy per nucleon, the pressure and 
the chemical potential, have jump discontinuities. The second order partial derivatives of the thermodynamical 
potential, i.e. the heat capacity per nucleon, the compressibility and the susceptibility per nucleon, have the 
infinite maximum or minimum at the points of phase transition. The phase diagrams T — v and T — p are 
represented by the coexistence lines; however, the phase diagrams T — s, T — p, T — /x, v — s, v — p and v — fx 
are depicted by the coexistence areas. The excitation energy per nucleon in the caloric curve is discontinuous in 
the canonical ensemble at constant value of temperature and fixed values of the specific volume. The pressure 
in the isotherms p — p and the chemical potential in the dependence p — p are discontinuous at constat value of 
the specific volume and fixed values of the temperature. 

We have shown that the first order phase transitions can be classified by the thermodynamic potentials of the 
statistical ensembles. We have also found that the first order phase transition associated with the Helmholtz 
free energy in the canonical ensemble requires the discontinuity of not only the extensive variables but also the 
intensive one. The Gibbs phase rule for this phase transition was formulated. It differs from the usual Gibbs 
phase rule of the first order phase transition of the liquid-gas type associated with the Gibbs free energy. The 
first order phase transition associated with the Helmholtz free energy is a consequence of thermal equilibrium 
between the system of nuclei and the environment which can exchange only heat with neither particle exchange 
nor work. 
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